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Abstract 

Carbon nanotube mechanical resonators have attracted considerable interest be¬ 
cause of their small mass, the high quality of their surface, and the pristine electronic 
states they host mm- However, their small dimensions result in fragile vibrational 
states that are difficult to measure. Here we observe quality factors Q as high as 5 x 10 6 
in ultra-clean nanotube resonators at a cryostat temperature of 30 mK, where we de¬ 
fine Q as the ratio of the resonant frequency over the linewidth. Measuring such high 
quality factors requires both employing an ultra-low noise method to detect minuscule 
vibrations rapidly, and carefully reducing the noise of the electrostatic environment. 
We observe that the measured quality factors fluctuate because of fluctuations of the 
resonant frequency. The quality factors we measure are record high; they are com¬ 
parable to the highest Q reported in mechanical resonators of much larger size [5J, [6|, 
a remarkable result considering that reducing the size of resonators is usually con¬ 
comitant with decreasing quality factors. The combination of ultra-low size and very 
large Q offers new opportunities for ultra-sensitive detection schemes and quantum 
optomechanical experiments. 
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In recent years, endeavours to boost quality factors in nano and micromechanical res¬ 
onators have been stimulated by the need to develop innovative approaches to sensing [7j, 
signal processing jS] and quantum physics [9]. Strategies to enhance quality factors have 
proceeded along three main routes. Firstly, the quality of the host material has been im¬ 
proved. To this end, new materials have been employed, such as high tensile stress silicon 
nitride membranes |5j HOJ and single crystal diamond films !H3- In addition, surface fric¬ 
tion has been lowered by optimizing fabrication processes and reducing contamination [12j. 
Secondly, schemes to isolate the resonator from its surrounding environment have been de¬ 
veloped, based on new resonator layouts H31 HU, and on optical trapping of thin membranes 
and levitated particles na m\- Thirdly, and most straightforwardly, Q-factors have been 
improved by operating resonators at cryogenic temperatures m 

Schemes to enhance Q-factors in nanotube resonators have focused on reducing contami¬ 
nation by growing ultra-clean nanotubes, and cooling resonators down to millikelvin temper¬ 
atures |2]. Even though Q-factors, measured from the linewidth of driven resonances, have 
been improved up to ~ 1.5 x 10 5 , they are still much lower than values routinely obtained 
with larger resonators fabricated from bulk materials using top-down techniques (6j. This 
result is somewhat disappointing, since the high crystallinity of nanotubes and their lack 
of dangling bonds at the surface are expected to minimize surface friction that limits the 
Q-factor in some nanomechanical systems [ IS] . 

Here we find that the actual values of the Q-factors can be significantly higher than 
hitherto appreciated, but that revealing these values requires to perfect the measurement 
technique. Namely, the dynamics of the nanotube has to be captured in a regime of van¬ 
ishingly small displacement in order to minimize nonlinear effects. In addition, noise from 
the electrostatic environment has to be reduced as much as possible. Indeed, owing to the 
ultra-small mass of nanotubes, this electrostatic noise affects the frequency of the nanotube 
resonator enormously and broadens the mechanical linewidth. All these experimental re¬ 
quirements make it challenging to unmask the instrinsic Q-factor. Progress in measuring 
ring-down with nanotubes has been made |19j . but it has not revealed higher Q-factor values; 
one reason might be that the large displacements in those measurements lead to sizeable 
nonlinearities. 

The geometry of our device, along with its characterization, are shown in Fig. 1. The 
nanotube is contacted by source (S) and drain (D) electrodes, and is suspended over a trench; 
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at the bottom of the trench is a local gate electrode (Fig. la). In this scheme, the ultra- 
clean nanotube is grown by chemical vapour deposition in the last step of the fabrication 
process of the resonator, making it free of fabrication residues. The device is cooled to 
30 mK, at which temperature all the data presented here are taken. We identify the lowest 
lying flexural mode of the resonator using the frequency-modulated (FM) mixing technique 
EH EDI, and study its dependence on gate voltage \V® ( | (Fig. lb). Figure lc shows resonance 
lineshapes in response to oscillating electrostatic forces, measured with the FM technique; 
they yield Q = 1700. The voltages used to produce these driving forces are kept low so that 
the resonance lineshapes are just above the noise floor. 

In the following, we employ an ultra-sensitive detection method described in Ref. [2] 
to capture the tiny amplitude of the thermal vibrations at 30 mK. In contrast to the FM 
technique, the resonator is not driven by an oscillating force. Displacement fluctuations are 
transduced into current fluctuations, of which power spectra Sj are measured (see Methods). 
We operate the resonator within the electron Fabry-Perot regime (Fig. Id), where the effect 
of electron transport on mechanical vibrations is less pronounced than in the Coulomb 
blockade regime (Fig. le) [2] 111 [23]. A typical resonance lincshape, obtained by averaging 
power spectra for ~ 512 s, is shown in the upper panel of Fig. 2a. The corresponding Q- 
factor is Q ~ 4 x 10 5 for a resonant frequency /o = 55.6 MHz. Remarkably, such a Q-factor 
is 200 times higher than that measured with the FM mixing technique (Fig. lc). 

We find that both the Q-factor and / 0 fluctuate in time. To show this, spectra obtained 
over a measurement time r = 3.2 s are acquired successively, using either a commercial source 
by Keithley or a simple lead battery to supply V g DC (see colour plots of Si as a function 
of frequency and time in Figs. 2a,b). Frequency fluctuations with the Keithley source are 
larger than with the lead battery. The Q-factor is fluctuating as well; its distribution is broad 
and asymmetric (Figs. 2c,d). On average, Q-factors are higher with the lead battery than 
with the Keithley source. While fluctuations in / 0 have been observed in liigh-Q mechanical 
resonators [2H - 126] . fluctuations in the measured Q-factor have not been discussed thus far. 
We show below that both fluctuations are closely related. 

The averaged Q-factor decreases with the measurement time (Fig. 3a). Once the mea¬ 
surement time is set, the Q-factor does not reveal any marked dependence on the amplitude 
V^ c of the oscillating source-drain bias used to read out current fluctuations (Fig. 3b). 
This shows that V^ c does not affect the measured motion of the resonator. For example, 
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if the resonator were to heat up due to Joule heating, the Q-factor would decrease as V^° 
increases; if the resonator were to sclf-oscillate, the Q-faetor would vary as V s ^ c increases. 
Neither the Q-factor nor / 0 changes as V g DC is swept through a conductance oscillation (see 
Supplementary Section VI). This confirms that the effect of electron transport on the vibra¬ 
tions is weak. This is crucial to our experiment, since Coulomb blockade would reduce the 
Q-factor mm- 

The observed frequency fluctuations are associated with the electrostatic noise of the 
environment. In the case of the Keithley source, frequency fluctuations are traced back to 
the voltage fluctuations of the source itself as a result of the strong gate voltage dependence 
of / 0 . To show this, we transform voltage fluctuations into Allan deviation a a using the /o- 
to-V g DC conversion factor measured in Fig. lb, ~ 1.4 x 10' Hz/V (see Methods). In Fig. 3c, 
we find that the measured Allan deviation of the resonant frequency of the nanotube (green 
hollow circles) coincides with the Allan deviation expected from the gate voltage fluctuations 
generated by the source (blue trace). The frequency noise is of l//-type (pink noise), since 
a a is essentially constant as a function of measurement time r. By contrast, the lead battery 
is stable enough that it does not significantly affect /o- The Allan deviation expected from 
the voltage fluctuations of the lead battery is smaller than a a(t) measured with the nanotube 
resonator (Fig. 3d). The origin of the observed frequency fluctuations with the lead battery 
is not clear. Since a a remains constant as a function of r (1// noise), / 0 fluctuations might 
be related to two-level systems such as charge fluctuators in the gate dielectric. 

The Q-factor is affected by the electrostatic noise of the environment as well. This can be 
inferred from a comparison between Fig. 2a and Fig. 2b. Larger fluctuations of the resonant 
frequency (lower panels) result in a larger resonance linewidth (upper panels), hence in a 
lower Q-factor. 

The fluctuations of the quality factor, the asymmetry of its distribution, and its depen¬ 
dence on measurement time are all attributed to the fluctuations of the resonant frequency. 
To show this, we separate this frequency noise into a slow part and a fast part on the scale of 
the ring-down time (see Methods and Supplementary Section X). For a measurement time r 
short compared to the characteristic time scale of the slow frequency noise t^ [ ow , resonances 
have a Lorentzian lineshape whose width 2r is constant while its resonant frequency fluc¬ 
tuates from one measurement to the next. For r t c s1ow , resonances are broadened by the 
time average of many such Lorentzians, resulting in non-Lorentzian lineshapes (Fig. 3e). For 
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r ~ Tg low , linewidths are larger than 2r, and fluctuate from one measurement to the next. 
Upon increasing r, linewidths tend to increase because the resonator has time to explore a 
larger frequency range. The distribution of linewidths is asymmetric since linewidths have 
a fixed lower bound given by 2T (see Supplementary Section X). 

Figures 2e,f display two high Q resonance lineshapes, obtained with both gate voltage 
sources, and demonstrating Q = 3.5 x 10 6 (Keithley source) and Q = 4.8 x 10 6 (lead 
battery). We verify that the areas of the resonances in Sj spectra in Figs. 2e,f are equal 
(within 10%) to the areas of the averaged resonances in Figs. 2a,b. This confirms that these 
sharp resonances do capture the dynamics of the nanotube, in spite of the relatively low 
number of measurement points dictated by the short measurement time. A slightly lower 
quality factor is obtained in an additional device, also operated in the Fabry-Perot regime 
(see Supplementary Section IX), confirming the robustness of high Q resonances. These Q- 
factors are comparable to the highest ones measured in large micromechanical resonators j3j. 
Yet, they may still be limited by the measurement, and may not have reached the intrinsic 
Q-factor of nanotube resonators, defined as /o multiplied by the small-amplitude ring-down 
time (the ring-down time is not known, but it is smaller than r; see Supplementary Section 
X). 

The large quality factors observed here are associated with the small amplitude of the 
nanotube vibrations. Using the equipartition theorem, we obtain that the mode temperature 
is T = 44 ± 10 mK, which corresponds to a phonon population n = ksT/hf 0 = 16 ± 4 with 
ks the Boltzmann constant and h the Planck constant (see Supplementary Section V). In 
this estimate, we use the areas of the resonances in Fig. 2, which are all equal within 10%. 
From the temperature, we calculate that the variance of displacement is ~ (35 pm) 2 . 

Larger displacements translate into lower Q-factors. This is illustrated in Fig. 3f where 
a white voltage noise is applied to the gate electrode in order to enhance the displacement. 
Possible origins of this behaviour could be associated to nonlinear damping forces, which 
result in a mechanical lincwidth that depends on the amplitude of motion, and to the spectral 
broadening of the resonance, which is induced by the combination of nonlinear conservative 
forces and displacement fluctuations [21 EM9]. Note that the applied voltage noise also 
induces fluctuations of / 0 , but their contribution to the measured resonance linewidth is 
negligible (see Supplementary Section VII). The temperature dependence of the Q-factor 
could not be measured in our current measurement setup. 
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The giant quality factors and the associated weak fluctuations of / 0 hold promise for 
diverse sensing experiments. The limit to force sensing is ultimately set by the force noise 
Sf = 87tM/c b T/ 0 /Q [2], Using M = 4.4 x 10” 21 kg, Q = 4.8 x 10 6 , and T ~ 44 niK, we 
obtain y/Sf ~ 1CU 21 N/\/Hz, which is lower than what has been achieved with mechanical 
resonators thus far [2j. The sensitivity of mass sensing and force gradient detection is given 
by the Allan deviation of frequency fluctuations (Sfo ~ 8 Hz for r ~ 6 s) in Fig. 3d obtained 
with the lead battery. The latter translates into a mass resolution SM = 2 M x 5/o//o — 
1CU 2 ' kg, which is as good as the best estimates on record [Tj. We extract a force gradient 
resolution of 8n 2 M f 0 Sf 0 ~ 1CU 10 N/m; this compares favorably with the best values reported 
thus far [30] . 

Nanotube resonators are relevant candidates for exploring the quantum regime of nanome¬ 
chanics. Their low mass leads to large zero-point motion, y/h/(AnfoM) ~ 6 pm, and al¬ 
lows to enhance the coupling to other quantum systems. High Q-factors are necessary for 
the quantum manipulation of mechanical states. In an experiment whereby the resonator 
is cooled to the ground state with passive feedback, the lifetime of the ground state is 
r p h = hQ/ksT ~ 10 -3 s (for ksT > hf 0 ) using T ~ 44 mK and assuming that Q = 4.5 x 10 6 
is related to the ring-down time. Our nanotube resonator demonstrates Qfo = 3 x 10 11 Hz, 
which compares well with the highest values measured thus far in nano and micromechanical 
resonators mm- 

Overall, our work shows that the Q-faetor of ultra-clean nanotube resonators can reach 
very high values, provided that the resonators are cooled to low enough temperature for 
nonlinear effects to be negligible, that slow frequency noise is reduced, that measurement 
times are short, and that the effect of the coupling between vibrations and Coulomb blockade 
is suppressed. The unique combination of ultra-low mass and giant quality factors offers new 
opportunities for ultra-sensitive detection schemes and optomechanical experiments in the 
quantum regime. 

Methods 

Detecting thermal vibrations. We employ an ultra-sensitive detection method that 
allows us to capture the tiny amplitude of the thermal vibrations in a dilution refrigerator 
cooled to 30 mK. Displacement fluctuations 8z in the direction z normal to the gate are 
transduced into conductance fluctuations SG = where C g is the capacitance 
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between the nanotube and the gate, and C g = dC g /dz (for the estimation of C g , see Supple¬ 
mentary Section II). Large transduced signals require a large transconductance dG/dV g . The 
largest transconductance is obtained in the Coulomb blockade regime, for V^ c > 0 (Fig. le), 
as is usually observed in ultra-clean suspended nanotubes [ 32 ]. However, Coulomb blockade 
reduces Q-factors, since it amplifies the coupling between electron transport and mechanical 
vibrations [ 21 - 01123 ] . For this reason, we operate the resonator within the Fabry-Perot regime, 
realized for V g DC < 0, even though dG/dV g is lower (Fig. Id). In this regime, oscillations 
in conductance originate from quantum interferences of electronic waves, and the effect of 
the coupling between vibrations and Coulomb blockade is weaker. The tiny conductance 
fluctuations 5G are parametrically down-converted into low frequency current fluctuations 
51 by applying a small oscillating voltage of amplitude V^ c across source and drain, at a 
frequency f s d shifted by ~ 10 kHz from the mechanical resonant frequency /q. Power spectra 
Sj of fluctuations 51 are measured at frequencies \f s d — fo\ ~ 10 kHz with a cross-correlation 
technique using a vector signal analyzer pj. 

Allan deviation. We calculate the Allan deviation a a as 


4(r) 


1 ( (fi+ i)t — ( fi)r \ 2 

2(A - 1) ^ V </o> ) 


( 1 ) 


where (fi+i) T and ( fi) r are two subsequent measurements of /o averaged over the integration 
time r, N is the number of averaged frequency measurements, and (/o) is the average of /o 
over the whole measurement [33]. We define / 0 as the frequency for which Sj is largest. 

Power spectrum of displacement and frequency noise. We separate the frequency 
noise into a slow part and a fast part on the scale of the ring-down time. The power spectrum 
of displacement S q (which is proportional to Sj ) as a function of angular frequency u = 2irf 
reads: 


S q (u) = [ dt- -—-, 

Muj 0 r J 0 T 2 + [u — uj 0 — £si ow {t )} 2 


( 2 ) 


where ks is the Boltzmann constant, T is the mode temperature, t is the time, and M is 
the modal mass. T is the sum of the reciprocal ring-down time and of the broadening due 
to fast frequency noise, Cjq is the resonant angular frequency renormalized by fast frequency 
noise, and C, s i ow is the slow part of frequency noise (see Supplementary Section X). The area 
Jo°° S q (cu)duj = nksT /Mojq is given by the equipartition theorem, hence it is independent of 
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Figure 1: Carbon nanotube mechanical resonator, (a) Schematic of the device. The nanotube 
is contacted by source (S) and drain (D) electrode, and is suspended over a gate electrode (G). The 
trench has a width of 1.8 gm and a depth of ~ 350 nm. (b) Resonant frequency /o as a function 
of gate voltage | V ; ^ c | for the lowest lying mechanical mode, (c) Mixing current I m i x as a function 
of drive frequency A/ (measured from /o = 44.1 MHz), using the FM technique. V^ c = —2.2 V. 
Two different driving voltage amplitudes are used: 13 (filled circles) and 8 fi\ (hollow circles). 
We use a lock-in amplifier with a time constant of 300 ms. In the FM technique, the Q-factor is 
extracted from the width of the foot of the resonance IS ESDI, yielding Q ~ 1700. (d-e) Conductance 
G as a function of Vg , in the Fabry-Perot regime (d), and in the Coulomb blockade regime (e). 
The full G{y^ c ) trace is shown in Supplementary Section I. 
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Figure 2: Fluctuations of the Q-factor and of the resonant frequency, (a, bottom) Power 
spectra of current fluctuations Sj(f — (/o)) acquired successively in time, using a commercial DC 
voltage source (Keithley) to supply the gate voltage {V® c = —3.037 V). For each of these spectra, 
the measurement time is 3.2 s. (a, top) Resonance obtained by averaging all the spectra, (b) Same 
experiment as in (a), this time using a lead battery to supply the gate voltage. Scale bar: blue, 
power spectral density Si = 0.15 x 10 -24 A 2 /Hz; red, Si = 2.50 x 10~ 24 A 2 /Hz. (c, d) Histograms 
of Q-factor using the Keithley source (c) and the lead battery (d), constructed from 3000 power 
spectra acquired with the same settings as in (a, b). (e, f) Examples of high Q resonances obtained 
with the Keithley source (e) and the lead battery (f). Red curves are fit to Lorentzian functions. 
For all panels, (/o) = 55.6 MHz is the time-averaged resonant frequency. 
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Figure 3: Characterization of the Q-factor and of the fluctuations in /q. (a) 1/Q as a 

function of measurement time r, using the Keithley source. For each measurement time, 300 power 
spectra are acquired and a histogram of 1/Q is built. Circles denote 1/Q at the maximum of the 
histograms, and bars represent the full width at half maximum of the histograms. The resonant 
frequency is /o = 55.6 MHz for panels (a-d). (b) 1/Q as a function of V//f' for t = 16 s. (c, d) Allan 
deviation of /o as a function of r with the Keithley source (c) and with the lead battery (d). The 
green circles are for the measured resonant frequencies, and the blue traces are for the resonant 
frequencies estimated from gate voltage fluctuations, (e) Schematic of resonance lineshapes for 
short (left) and long (right) measurement times r compared to the characteristic time of the slow 
frequency noise t/ Iow . ii, 2,3 are three successive instants, (f) 1/Q as a function of added gate 
voltage noise at V^ c = —2.2 V (/o = 44.1 MHz) for r = 16 s. The gate voltage noise is applied 
using the Johnson-Nyquist noise of a 50 Ohm resistor at room temperature amplified by different 
gains. It creates a white electrostatic force noise between the nanotube and the gate (top axis). 
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I. CONDUCTANCE AS A FUNCTION OF GATE VOLTAGE 


Fig. 


The conductance G as a function of gate voltage V® c at 30 mK is shown in Supplementary 
The regime for positive V g corresponds to Coulomb blockade (Supplementary 


SI 


Fig. SI 3, middle and right-hand-side panels). Indeed, the length of the blockaded region, 
estimated to be ~ 1.5 pm from the separation between conductance oscillations, is close to 
the width of the trench (see next section). For negative V g DC , charge transport is in the 


so-called Fabry-Perot regime (Supplementary Fig. SI 3, left-hand-side panel). The period in 
V g DC of the oscillations of G is a few times larger than the period of oscillations for positive 
V g DC , as usually observed for suspended nanotubes. In addition, the conductance is close to 
2e 2 //i, indicating that the electronic transmission of each contact is high. 


The measurement shown in Supplementary Fig. [ST] is typical of ultra-clean nanotubes 
[1]. The regions of the nanotube near the metal electrodes remain p-doped, due to the work 
function of the electrodes, whereas the suspended part of the nanotube can be doped with 
electrons or holes using V g DC . For positive V g DC , p — n junctions are formed near the metal 
electrodes, resulting in a Coulomb blockaded region along the suspended nanotube. For 
negative V ® c , the nanotube is p-doped along the whole tube and no tunnel barriers are 
formed. Near V g = 0, the conductance is zero because the Fermi energy lies inside the 
energy band gap along the entire suspended part of the nanotube. 



Figure SI: (a) Conductance G as a function of gate voltage V g DC • Regions within colored rectangles 
are blown up in (b). 
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II. ESTIMATING C' g 


We estimate the capacitance C g between the nanotube and the gate from the separation 
AVg = 13 ± 1 mV between two conductance peaks in the Coulomb blockade regime 
(Supplementary Fig. SI 3, central and right-hand-side panels): C g = e/AV^ C = 1.2 ± 0.1 x 
10 -1 ' F. This value is close to the model capacitance between a cylinder of length L and 
radius r, and a plane a distance h away: 


C a 


2Tre 0 L 

In ’ 


where eo is the vacuum permittivity, h is the separation between the nanotube and the gate, 
and z ^ h is a small displacement of the cylinder in the direction normal to the gate. The 
derivative dC g /dz reads 




dC n 


c n 


= (5.2 ± 1.9) x 10~ 12 F/m, 


dz J z=0 h\n(2h/r) 

using h = 350 ± 50 nm and r = 1 ± 0.5 nm (the typical radius of nanotubes obtained with 
our chemical vapour deposition technique). 


III. ESTIMATING THE EFFECTIVE MASS 

We study the fundamental flexural mode of the nanotube resonator (see Fig. lb of the 
main text). We use a simple model whereby the nanotube is straight and is perpendicular 
to the trench. The effective mass M of the mode is related to the mass of the nanotube 

Vf/VT 8S 

1 f L 

M = M nt — / [(f)(x)] 2 dx , 

L Jo 

where L is the length of the nanotube and (f>(x) is the shape of the mode, which is normalized 
so that max[0(x)] = 1. Given that mechanical tension is induced in the nanotube by the 
contacts and by the gate voltage, we assume that the modal shape is (p(x) = sin(7r x/L). The 
latter is the simplest approximation for the shape of a beam under tension. We emphasize 
that the expression for M takes into account the shape of the vibrational mode; all other 
quantities are measured with respect to the amplitude of this mode. To estimate M NT , we 
also assume that the length of the nanotube is equal to the trench width (L = 1.8 ± 0.2 /jm) 
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and that r = 1 nm as in Section II. Therefore 


j = (4.4 ±2.4) x KT 21 kg, 

where Me = 2 x 10~ 26 kg is the mass of a carbon atom and A = 5.2 x 1CF 20 m 2 is the surface 
area of a hexagon in the honeycomb lattice of graphene. 


M — - 
2 


2 Me x 


2nr x 


A 


IV. ORIENTATION OF THE MODE WITH RESPECT TO THE GATE ELEC¬ 
TRODE 


The orientation of the vibrations of the eigenmode depends on the static curvature of 
the suspended nanotube, which builds in during the fabrication of the resonator. Since the 
static curvature cannot be controlled and cannot be accurately measured, the eigenmode 
can vibrate in any direction. For a real nanotube, the displacement along the nanotube can 
be quite complicated. However, we will use a simplified model where at least for the lowest 
mode the displacement is in one plane. We will characterize this plane by the angle 6 it 


makes with the plane parallel to the surface of the gate electrode (Supplementary Fig. S2i). 

We assume that current fluctuations 51 at the drain electrode are proportional to motional 
fluctuations Sz along the direction z normal to the gate electrode. Then, the current at the 
frequency close to the difference between the mode eigenfrequency and the frequency of the 
source-drain voltage is 

1 rin __ . CV 

(SI) 




where dG/dV g is the transcondutance, V^ c is the static gate voltage, and V^j c is the am¬ 
plitude of the oscillating source-drain voltage. We assume that the mode is polarized along 


q, as shown in Supplementary Fig. S2i; hence, 6z is the projection along z of the motional 
fluctuation 5q along q: 

8z = 5qsm6 , (S2) 


where 9 is the angle between y and q. To estimate 0, we subject the nanotube to a weak 
electrostatic force by applying a small oscillating voltage with amplitude 8V^ C to the gate. 
The component of the force along z is 8F Z = CgV^ 0 5V^ C . The center of mass of the 
nanotube experiences a force 8F q = 6F z sin9, the projection of 5F Z along q. On resonance, 
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a 


b 




Figure S2: (a) Orientation of the displacement for a simple model of a straight nanotube parallel to 
the gate electrode, x and y are parallel to the gate, and z is normal to the gate. The displacement 
(that is, the mode polarization) is along q. The force SF Z along z acts on the motion of the 
nanotube through the force 5F q , its projection onto q. The projection of displacement 5q onto 
z is 5z. (b) Standard deviation of current yj(SI 2 ) induced by a small oscillating gate voltage of 
amplitude SV AC applied on resonance. SV AC produces the force 5F Z in (a). 


the displacement Sq induced by SF q reads 


5q = 


-0—5F 
Muf q ' 


(S3) 


where Q is the quality factor of the resonance, M the effective mass, and the resonant 
frequency. In turn, the displacement Sq induces current oscillations of variance (SI 2 ) at the 
drain: 


VW) = P sin 2 0-^C'V DC $ V q AC = pSV AC . 


Mu 2 9 ' 9 


(S4) 


Supplementary Fig. S2 i shows \J(SI 2 ) as a function of SV AC , measured with the technique 
described in the main text and in Ref. [2]. The dependence of a/ (SI 2 ) on SV g AC is linear 
and extrapolates to zero in the limit of zero SV AC , as expected from Eq (S4). We obtain 
6 = 60° ± 25° from the slope in Supplementary Fig. S2 3, using \V g Di '\ = 2.2 V, V S AC = 
3.17 x KT 4 V, C g = (1.2 ± 0.1) x 10" 17 F, C' g = (5.2 ± 1.9) x 10~ 12 F/m, M = (4.4 ± 
2.4) x 10~ 21 kg, wo/(27t) = 44.1 x 10 6 Hz, Q = 5 x 10 5 , dG/dV a = 6.4 x 10“ 4 S/V, and 
p= (1.3 ±0.2) x 10" 3 A/V. 

The estimation of 9 comes with a large uncertainty that mostly originates from the 
uncertainties in the estimations of C' g and M. These uncertainties, however, have little 
impact on the estimation of the modal temperature T, which we present is Section V. 
Indeed, we show that T is proportional to M/(C' g sin 9) 2 , a quantity that we can estimate 
rather precisely from the measurement shown in Supplementary Fig. 2b. For this, we rewrite 
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Eq. (S4) as: 


M 


(C' smO) 2 p \2 dV 6 


— 1 / I dG n/ -pc\2-[/-AC -*■ Q 

_ ) V SCL CgUjl 


(S5) 


Using the values given above, we calculate the right-hand-side of Eq. (S5) and obtain 

M = 205 ± 40 kg ■ m 2 /F 2 . 


(C' sin 9) 2 


We use this value to estimate T is Section V. 


V. MODE TEMPERATURE 


Using Eqs. (SI) and (S2), the equipartition theorem yields 

1 I r°° 

k B T = Mul(5q 2 ) = Mu 2 0 


Si(u)duj , 


(S6) 


°{P sm9) 2 2n J 

where ( Sq 2 ) is the variance of displacement along q, and Sj(lo) is the two-sided power spectral 
density of current fluctuations 61. To make contact with the experiment, we also express 
k B T in terms of Si(2nf), where f — is the natural frequency: 


k B T = Moj 2 


(£sin 9) 2 J 0 


2S T (2nf)df, 


(S7) 


where we used 


— S'/(a;)da; = / Sj{2^f)df = 2 / S z (2nf)df. 

27T ,/n 


In Eq. (ST), Si{2nf) is defined as 


r-r/2 


Si(2nf) — lim / (<5/(f)<5/(0)) exp[-i(27r/ + o; sd - o; 0 )f]df, 


(S8) 


'-t/2 


where r is the measurement time, uj s d is the off-resonance angular frequency of the source- 
drain voltage, and ujq is the resonant angular frequency. Experimentally, we measure the 
single-sided power spectral density S'| xp (27r/) = 2S/(27 t/). 

We estimate the mode temperature T from the power spectra in Figs. 2e,f of the main text. 
We obtain T ~ 44 ± 10 mK, using the following parameters: \V g DC \ = 3.037 V, dG/dV g = 
8.5 x 10" 4 S/V, Vff = 4 x 10- 4 V, C g = (1.2 ± 0.1) x 10~ 17 F, o; 0 /(27r) = 55.6 x 10 6 Hz, 
and / 0 °° *S'/ XP (2vr/)d/ = 4.5 x 10~ 23 A 2 . We also use M/(C' g sin 9) 2 = 205 ± 40 kg ■ m 2 /F 2 


obtained in Section IV. 

We were not able to carry out a temperature dependence of ( Sz 2 ). Connecting thermom¬ 
etry lines to the cryostat would generate electrical noise in our measurement, masking the 
mechanical resonance in the spectra. 
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VI. DEPENDENCES OF THE RESONANT FREQUENCY AND THE Q-FACTOR 
ON DC GATE VOLTAGE 


Supplementary Fig. S3i shows the resonant frequency as a function of V® c obtained by 
measuring the mixing current I m i x with the FM technique [3] as a function of V g DC and 
drive frequency /. The intensity of I mix varies with V g DC because it is proportional to the 
transconductance. The resonant frequency / 0 is linear in V g DC . Supplementary Fig. S3 
shows the inverse of the Q-factor as V g is stepped from the bottom of an oscillation 
I DC (V DC ) to the top of this oscillation. No variation of 1 /Q is 


m 


seen. 


The dependences of /o and 1/Q on V g are different from what is observed in nanotube 
resonators operated in the Coulomb blockade regime. There, the resonator experiences a 
reduction of /o near the top of an oscillation in I DC as a function of V® c '. In addition, 
1/Q is higher near the top of the oscillation in I DC because electron tunneling events are 
accompanied by enhanced dissipation. The absence of such behaviors demonstrates that the 
effect of the coupling between vibrations and Coulomb blockade is weak in our experiment. 


VII. EFFECT OF ELECTROSTATIC WHITE NOISE ON THE SPECTRUM OF 
THE RESONATOR. 

Figure 3f of the main text displays 1/Q in the presence of gate voltage noise of various 
intensities. The gate voltage noise is applied using the Johnson-Nyquist noise of a 50 Ohm 
resistor at room temperature amplified by different gains (by varying the number of ampli¬ 
fiers). We measure the power spectral density of the amplified Johnson-Nyquist noise with 
a signal analyzer in a separate experiment. We verify that this amplified noise is white (the 
power spectral density is constant) between 1 MHz and 200 MHz, a frequency range that en¬ 
compasses the resonant frequency of our resonator (40-60 MHz). Noise at frequencies below 
1 MHz is cut off by a high pass filter at the sample stage. The power spectral densities we 
measure are consistent with the gain (20 dB) and the noise figure (1 dB) of each amplifier 
used in combination with the total attenuation along our radio frequency lines. 

The applied Johnson-Nyquist noise creates a random electrostatic force between the nan¬ 
otube and the gate (top axis of Fig. 3f). To the lowest order in the nanotube displacement, 
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Figure S3: (a) Electro-mechanical mixing current measured with the FM technique, as a function 
of V g DC and drive frequency /. Scale bar from blue, 1 pA, to orange, 140 pA. (b) Source-drain 
DC current I DC as a function of V g , showing Fabry-Perrot oscillations, (c) 1/Q as a function of 
time averaged V^ C ■ Measurement time is 6.3 s. The error bars correspond to the distribution of 
Q-factors for 20 power spectra. 

this force is additive, which means that it is independent of the displacement. The nanotube 
resonator responds primarily to the frequency components of this force, which lie within the 
band centered at the resonant frequency, with the typical bandwidth given by the mechan¬ 
ical linewidth. The power spectrum of the random force is flat in this region, and therefore 
the corresponding noise is white. Its intensity S e f is related to the power spectral density 
of the Johnson-Nyquist noise Sjn as S e f = {C' g V^ c s\i\9) 2 Sjn (see Section IV). We use 
C' g = 5.2 x 10- 12 F/m, V g DC = 2.2 V, and 9 = 60°. 

The effect of the additive Johnson-Nyquist noise is fully analogous to the effect of the 
thermal noise that comes along with the resonator dissipation. It leads to the increase of 
the intensity of the peak in the resonator power spectrum and to spectral broadening of 
this peak via nonlinear damping and via the dependence of the oscillator frequency on the 
vibration amplitude. 

The low-frequency components of the Johnson-Nyquist noise could directly lead to fluc- 
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tuations of the resonator frequency, which would be similar to the low-frequency fluctuations 
induced by the voltage noise of the DC voltage source. They would come from the depen¬ 
dence of the resonator frequency on the gate voltage. Such noise is multiplicative, since the 
corresponding force is proportional to the resonator displacement, and this is why it acts 
as a frequency shift. The resonator responds to such noise if its frequencies are low, within 
the typical bandwidth given by the mechanical linewidth. Higher-frequency components are 
averaged out. For the Johnson-Nyquist noise that we are using the effect of this noise is 
very weak, because the initial noise from the 50 Ohm resistor is only weakly amplified below 
100 kHz, and is further attenuated by a high-pass filter below 10 kHz. We estimate that the 
corresponding gate voltage fluctuations induce a lineshape broadening of only a fraction of 
a Hertz, much smaller than the linewidths we measure in Fig. 3f. 


VIII. VOLTAGE FLUCTUATIONS OF THE DC VOLTAGE SOURCES 


The voltage fluctuations of both the Keithley source and of the lead battery are measured 
in a separate experiment. The DC source is connected to the DC gate voltage input port 
outside of the cryostat. A finite DC voltage close to the gate voltages used in the main 
text is applied. The voltage fluctuations at 300 K are measured at the sample stage after 
opening the cryostat. These fluctuations are recorded with a voltmeter as a function of time 


(Supplementary Fig. S4). The measured voltage fluctuations are larger than the intrinsic 


readout fluctuations of the voltmeter. These measurements are used to estimate the Allan 
deviation in Figs. 3c,d (blue traces) of the main text. 


IX. ADDITIONAL SAMPLE 


We investigated a second nanotube resonator device operated in the Fabry-Perot regime 
at a refrigerator temperature of 30 mK. An example of a high Q resonance for this device 

The parameters used for this measurement are: V® c = 


S5 


is shown in Supplementary Fig. 

—3.643 V, Vff' = 4 x 10 -4 V, and a measurement time of 4.8 s. The measurements were 
carried out with the Keithley source. 
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Figure S4: Gate voltage fluctuations measured at the sample stage, (a) Keithley source, (b) Lead 
battery. 



Figure S5: High Q resonance obtained with an additional device. The resonant frequency is 
f 0 ~ 62 x 10 6 Hz. 


X. POWER SPECTRUM FOR A FINITE MEASUREMENT TIME 


Below, power spectral densities are all two-sided. 
fluctuations for a finite measurement time r reads 

S q (u}) — — [ dt f dt\ q(t)q(ti)e luJ ^ tl ' > = -Re 
T Jo Jo T 


The power spectrum of displacement 



dtiq(t)q(ti)e 




(S9) 
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We make the change from fast oscillating variables q(t), q(t ) to slow complex oscillator 
amplitude u(t) 


q(t) = u(t)e l “ ot + c.c., 
q(t) = ioo 0 {ue iuJot - u*e~ iuJot ) . 


(S10) 


Plugging Eqs. (S10) into the equation of motion q + 2Tq + [oJq + 2uo£(t)]q = frift/M, and 
solving for u in the rotating wave approximation, we obtain 


u(t) — / dt' exp 


-T{t-t')+i / £{f)df 


Ut'). 


(Sll) 


In the above expression, £ is the frequency noise, f u = l/(2iMcu 0 )fT(t) exp(— icuot) is the 
thermal noise with correlator (f u (t)f u *(t')) = v/fj ^~ *0, an d M is th e m °dal mass. Using 


Eqs. (S10) and (Sll) along with Eq. (S9), we obtain 


SJu) — —Re dt dt\ u*(t)u{ti) exp [i(u — u>o)(t — U)] 

T ,/n ./n 


rt 


rti 


= -Re dt dt\ 

T I n ,/n 


/ dt ' / 

oo J — oo 

rt 1 


x exp 


Tk B T2 


Mojq t 


-i / mdt "+* / mdt'i 


dt[ exp [i{uj — uo)(t — t\) — T(t — t') — T(ti — ^)] 
Tk B T 


x 


Mluq 


S(t' -1[) 


ft i 


-Re dt d.t\ / dt' l exp [i(uj — uio)(t — t±) — T(t + t\ — 2t')\ 


x exp 
k B T 1 


Mwg t 


at")dt‘ 


-Re dt dt] exp 


i(co ~ u>o)(t — ti) — T(t — t\) — i / ^{t")dt h 


't\ 


(S12) 


In Eq. (S12), we average over thermal noise since its correlation time is short on the scale 
r- 1 <C r. We can now separate £( t ) into two parts, one slow and the other fast on the scale 
of the ring-down time r -1 : £(t) = £, s iow(t) + £,fast(t)- We assume that £ f ast is (^-correlated, 


which allows to simplify the power spectrum in Eq. (S12) as 


k B T 1. 


S •<"> = Mug 


-Re dt dt i exp 


r 


k B T 1 


dt 


Muq t J 0 f2 + [tu-Cu 0 -Z slow (t)}■ 2, 


i(u -u 0 ) - T (t - ti) - i£siow{t)(t - ti) 
|(U — CUq | < ^. ^0 j 


(S13) 
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where f is the “instantaneous” half-width of the spectrum and loq is the resonant angular 
frequency, both renormalized by fast frequency noise. The relation between T and T can be 


found in Supplementary Ref. pE]. Equation (S13) is Eq. (2) in Methods. The resonant peak 
corresponding to S q {uj) only exists for exp frj <C 1; the fact that resonance lineshapes 
remain close to Lorentzian indicates that our measurement times exceed 1/T. 

We emphasize that £,f a st(t) is 5-correlated on the “slow” time scale ~ T _1 , not on the fast 


scale ~ 0 Jq . It is seen from Eq. (S12) that the components of £(t) with frequencies much 
higher than T are averaged out and therefore can be disregarded, as we indicated earlier. The 
case where the frequency noise has significant intensity near 2ujq, so that it parametrically 
excites the resonator, requires a separate analysis, but we have no indications and no physical 
reasons to expect that strong frequency noise with frequencies ~ 2wo is present in our case. 

The separation of the frequency noise into parts that are slow and fast on the time scale 
1/T leaves out a comparatively narrow part of the noise spectrum. This part is averaged out 
when, as in our case, the duration of a measurement r 1/T. It is reasonable to expect that 
the contribution of this narrow range of the frequency noise spectrum is small. Additional 
information about the spectrum of the frequency noise can be obtained by studying the 
power spectrum of the resonator in the presence of periodic modulation j?j. 


The integral width of the spectrum can be obtained from Eq. (S13). The integral width 
is expressed as the ratio X = A/S™ x (cu), where A = J 0 °° S q (co)du is the spectrum area and 
5“ ax (w) is the resonance height in the spectrum. The area A = nksT /(Mwq) is independent 
of frequency noise or decay rate. For a Lorentzian spectrum (£, s i ow = 0) we have 

X" 1 = i = Q-?-, 

7lT vtw 0 

where Q is the quality factor. 

Next, we calculate X -1 , a good number to characterize Q for asymmetric and noisy 
resonances. For slow frequency noise of weak intensity, the expression for the reciprocal 
width takes on a “finite-time variance” form: 

[ — (— f dt £slow(t) 


X -1 ~ Wr | 1 - ^ 


7t r 


r 2 


(S14) 


where ^ dt £, s i ow (t) is the measured shift of the resonant frequency resulting from slow 
frequency noise. The term in square brackets in Eq. (S14) is always positive, resulting in 
the two following effects: 
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Keithley source 


Lead battery 



Figure S6: Histograms of the reciprocal integral width X 1 , using the Keithley source (a) and the 
lead battery (b) to bias the gate electrode. 

1. For a fixed measurement time r, multiple measurements of the integral width give an 
asymmetric distribution of X -1 ; 

2. As r increases, the peak value of the distribution of X -1 decreases, whereas the mean 
X increases. Once r exceeds the noise correlation time, X levels off to the value 

(i+ £<&«,(*)>)■ 


Supplementary Fig. S6 shows histograms of X -1 built from the spectra used to construct 
the histograms in Figs. 2c,d of the main text. If, as we assume, the modulation of the 
conductance is linear in the resonator displacement, the current-to-displacement scaling 
coefficient drops out from X. Therefore, the value of X can be found directly from the 
data on the current power spectra. Histograms of X -1 are asymmetric, in agreement with 
effect (1) predicted above. In addition, Fig. 3a of the main text shows that the peak value 
of the linewidth distribution tends to level off for long measurement times, in agreement 
with effect (2). (From that figure, we infer that the spectrum of the frequency noise has 
a cutoff around 1/20 Hz.) Slow fluctuations of the decay rate and/or the intensity of the 
thermal noise would not lead to the observed increase of the the peak value of the linewidth 
distribution with measurement time, but rather would cause X to decrease. 
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